flD-fll58  573  STOCHASTIC  REARRANGEMENT  INEQUALITIES(U)  ELORIM  STATE  T 
UNIV  TALLAHASSEE  DEPT  OF  STATISTICS  C  D^ABADIEL  ET  AL. 

SEP  82  FSU-STRTISTICS-M672  AFOSR-TR-85-0807 
UNCLASSIFIED  F49620-82-K-0007  F/G  20/1  NL 


AFOSRTR-  85-0  007 


STOCHASTIC  REARRANGEMENT  INEQUALITIES 


Catherine  D'Abadie^ 

Western  Electric  and  Florida  State  University 


Frank  Proschan^ 
Florida  State  Ikiiversity 


FSU  Statistics  Report  No.  H672 
AFOSR  Technical  Report  83>167 


September,  1983 


The  Florida  State  University 
Department  of  Statistics 
Tallahassee,  Florida  32306 


"tv.,  r.  t 


i  ^  Ld  i 


1  /' 

Research  sponsored  by  the  Air  Force  Office/of  Scientific  Research,  Air  Force 


Systems  Command,  USAF,  under  Grant  Number  AEOtR-82-K-0007.  The  U.S.  Government 
is  authorized  to  reproduce  and  distribute  reprints  for  Governmental  purposes 
notwithstanding  any  copyright  notation  thereon. 

Key  Words:  Rearrangement  inequalities,  stochastic  rearrangement  ineqxialities , 
decreasing  in  transposition,  ranking  problems,  optimal  assembly  of  systems, 
contamination  models,  hypothesis  testing,  partial  ordering,  total  positivity, 
positive  set  function,  arrangement  increasing. 


AMS  siibjeat  olaasifiaation:  60C05 


reieft* 


UNCLASSIFIED 


'If (CATION  Of  this  fACE 


AD-AIS0S73 


REPORT  DOCUMENTATION  PAGE 


1*  HEfOAT  SECUNITV  CLASSIf  (CATION 

UNCLASSIFIED 


ID.  ME5THICTIVE  MAAKINCS 


3»  SECUNITV  CtASSIf  (CATION  AUTHOAITV 


3b.  OECuASSIfJCATION/DOINNGMAOING  SCHEDULE 


3.  OiSTNISUTION/AVAILABILITT  Of  MEPOAT 

Approved  for  public  release;  distribution 
unlimited. 


A  PEAfOANHNfi  OACANIZATION  NEPOAT  NUMBEA(S) 

FSU  TR  No..  M672;  AFOSR  TR  No.  83-167. 


5.  MONITOAING  OAGANIZATION  AEPOAT  NUMBEA(Sl 

AFOSR-TR-  8  5-0007 


6*.  NAIME  Of  PEAfOAIMINC  OAGANIZATION 

Florida  State  University 


Of  f  ICE  SYMBOL 
(It  •ppUcmbl*  I 


T*.  NAME  Of  MONITOAING  OAGANIZATION 

Air  Force  Office  of  Scientific  Research 


6«.  AOOAESS  (City.  Slaf*  and  ZIP  Codtl 

Dept  of  Statistics 
Tallahassee  FL  32306 


7b.  ADOAESS  (City.  Stair  and  HP  Coda) 

Directorate  of  Mathematical  &  Information 
Sciences,  Bolling  AFB  DC  20332-6448 


la.  NAME  Of  f  UNDING/SPONSOAING 
OAGANIZATION 

AFOSR 


b.  Off  ICE  SYMBOL 
(It  appUeablai 

NM 


9.  PAOCUAEMENT  INSTAUMENT  lOENTIf  ICATION  NUMBEA 

F49620-82-K-0007 


Be  AOOAESS  (City.  Slate  and  ZIP  Coda) 


Bolling  AFB  DC  20332-6448 


10.  SOUACE  Of  f  UNDING  NOS. 


PAOGAAM 

PAOJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

61102F 

2304 

A5 

11.  title  (Ineluda  Sacurtty  Claaiifiaationi 

STOCHASTIC  REARRANGEMENT  INEQUALITIES 


13.  PEASONAL  AUTHOAIS) 

Catherine  D'Abadie  and  Frank  Proschan 


Il3*i  tvk*  Ob  PC*0«T 

13b.  TIME  COVEAEO 

1«.  DATE  Of  AEPOAT  (Yr.  Mo..  Day) 

15.  PAGE  COUNT  1 

1  Technical 

fAOM  TO 

SEP  84 

38  1 

IS.  SUPPLEMENTAAY  NOTATION 


17. 

COSATl  COOES 

18  SUBJ8CT  T8NMS  rConfifiM  on  pviwpvc  if  neceMory  uh4  identify  by  btoek  nvmBtr; 

flELO 

GNOUA 

SUB.  CA 

Rearrangement  ineo^alities;  stochastic  rearrangement 

inequalities;  decreasing  in  transposition;  ranking  problems; 

optimal  assembly  of  systems;  contamination  models;  CONTINUED 

19.  A9STMACT  (Continut  on  if  ntetmury  tdtntify  by  biocb  Humbert 

The  authors  develop  a  unified  theory  for  obtaining  stochastic  rearrangement  inequalities. 
The  authors  present  sample  applications  in  ranking  problems,  hypothesis  testing,  contamina¬ 
tion  models,  optimal  assembly  of  systems,  and  stochastic  versions  of  well  known  rearrange¬ 
ment  inequalities..,.^^ 


ITEM  //18,  SUBJECT  TERMS,  CONTINUED;  'hypothesis  testing;  partial  ordering;  total  positivity; 


positive  set  function;  arrangement  increasing. 

30.  OISTAIBUTION/AVAILABILITY  Of  ABSTAACT 

31.  ABSTAACT  SECUAITY  CLASSIf  ICATION  I 

UNCLASSIflEO/UNLIMITEO  E)  SAME  AS  APT.  □  OTIC  USSAS  □ 

UNCLASSIFIED 

33a.  NAME  Of  AESPONSiSLE  INDIVIDUAL 

33b  TELEPHONE  NUMBEA 
ilnaluda  Araa  Coda) 

33c  Off  ICE  symbol 

MAJ  Brian  W.  Woodruff 

(202)  767-  5027 

NM 

OD  FORM  1473,  83  APR  edition  of  i  jan  73  is  obsolete.  UNCLASSIFIED 


BECUAITY  CLASSIf  ICATION  Of  THIS  PA  \E 


Stochastic  Rearrangement  Inequalities 


by 

Catherine  D'Abadie  and  Frank  Proschan 


ABSTRACT 


We  develop  a  unified  theory  for  obtaining  stochastic  rearrangement 
inequalities.  We  present  san^le  applications  in  ranking  problems,  hypothesis 
testing,  contamination  models,  optimal  assembly  of  systems,  and  stochastic 
versions  of  well  known  rearrangement  inequalities. 
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1.  Introduction  and  SiinmaTy.  We  obtain  stochastic  versions  of  rearrangement 
inequalities.  Rearrangement  inequalities  compare  the  value  of  a  function  of 
vector  arguments  with  the  value  of  the  same  function  after  the  components  of 
the  vectors  have  been  rearranged. 

The  classical  example  of  a  rearrangement  inequality  involving  a  function 
of  2  vector  arguments  is  the  well-known  inequality  of  Hardy,  Littlewood,  Pc/lya 
(HLP)  (1952)  for  sums  of  products.  For  vectors  (aj.-.-.a^^)  and  fe*  (bj,...,bjj) 
of  positive  numbers,  HLP  show  that  if  a^i  ...  and  bj^i  ...  ab^^  (after  relabeling, 

say),  for  every  permutation  (Tr(l) , . . .  ,iT(n))  of  (l,...,n),  then 

n  n  b 

y  a.b.  i  y  a.b  ...  ^  I  a.b  .  , 

^,11  1  ir(i)  1  n-i+1 


i«l 


i*l 


i=l 


(1.2) 


r  St .  U  St .  S 

I  Vi  ‘  ‘ 

1«1  1»1  ^  1»1 


(1.1) 

hold. 

What  is  a  stochastic  version  of  (1.1)?  Under  what  conditions  will 

”  St.  2 

V  V  >  \  V  V 

t(i) 

hold,  given  random  vectors  X,  X’  To  answer  this  question,  we  need  the  following 
definition: 

^ii 

1.1  Definition.  For  vector  x  we  write  x  2  x'  if  i <  j ,  x.  s x .  and  x'  is  obtained 

I  ■  I  ■  ■  I  M  13^ 

from  X  by  interchanging  x.  and  x.  and  leaving  the  other  components  fixed. 

1.2  Proposition.  Let  X  and  Y  be  nonnegative  random  vectors  having  joint  density 
f(x;^).  Then  inequality  (1.2)  holds  for  }C  and  Y  if  for  all  pairs  (i,j),  isi<jsn, 
f  satisfies: 

^ij  ^ij 

(1.3)  f(x,;^)  ♦  f(x',x')  -  ^  ***'®’^®  ii  ^  X'  and  jr  2  y'. 

Many  well-known  multivariate  densities  satisfy  this  condition,  as  we  shall 

see. 


Since  the  work  of  HLP,  many  papers  on  rearrangement  inequalities  have  appeared. 


More  recent  are  London  (1970) ,  Mine  (1971) ,  and  Day  (1972) .  The  Marshall  and 


Olkin  (MO)  (1979)  book  presents  a  unified  approach  to  the  study  of  deterministic 
rearrangement  inequalities. 

In  this  paper  we  develop  a  unified  theory  for  obtaining  stochastic  versions 
of  rearrangement  inequalities.  (1.2)  becomes  a  special  case,  as  does  the  work 
of  many  earlier  authors.  Moreover,  we  obtain  stochastic  refinements  of  these 
inequalities  analogous  to  those  obtained  by  MO  in  the  deterministic  case. 

The  value  of  having  stochastic  versions  of  rearrangement  inequalities, 
apart  from  purely  mathematical  interest,  is  manifested  in  their  applicability 
in  a  surprisingly  large  number  of  statistical  and  reliability  contexts  such 
as,  e.g.,  ranking  problems,  hypothesis  testing,  contamination  models,  and 
optimal  assembly  of  systems. 

In  Section  2  we  present  some  definitions  and  concepts.  In  Section  3  we 
establish  preservation  properties  for  the  functions  of  interest  under  various 
statistical  and  mathematical  operations.  The  theory  of  stochastically  similarly 
ordered  random  vectors,  which  provides  a  unified  ^proach  to  stochastic  rearrange¬ 
ment  inequalities,  is  developed  in  Section  4.  In  Section  4  we  also  show  that 
many  well  known  multivariate  densities  govern  stochastically  similarly  ordered 
random  vectors.  In  Section  5  we  present  some  illustrative  applications  of  the 
theory  to  ranking  problems,  hypothesis  testing  problems,  and  contamination 
models. 

2.  Preliminaries.  In  this  section  we  introduce  definitions  and  preliminary 
results  used  in  subsequent  sections. 

Let  r”  denote  Euclidean  n-space.  Hfe  define  a  partial  ordering  of  r"  x 
which,  as  MO  (1979)  show,  unifies  the  study  of  deterministic  rearrangement 
inequalities.  To  define  this  partial  ordering  we  need  some  terminology  and 


notation. 


Let  denote  the  group  of  all  permutations  of  {1,  2 ,  . . . ,  n} .  An  element 
of  is  denoted  by  £=  (ir(l)  f  ....  ir(n)).  We  sometimes  identify  with  the 
subset  of  whose  elements  are  those  vectors  with  the  integer  components 
1,  2,  n.  Let  ir  and  it'  be  elements  of  S_.  We  say  that  it'  is  a  simple 

M  JJ  ^ 

transposition  of  ir  if  there  exist  positive  integers  1  s  i  <  j  s  n  such  that 

ir(i)  *  ir'(j)  =  Tr'Ci)  =  Tr(j)  and  ir(k)  *  ir'(k)  for  k  *  i,  j.  We  write  this  as 

^ii  t 

IT  >  tt".  For  IT,  it'  e  we  say  that  ir^  is  a  transposition  of  jr,  written  viv', 

if  w  s  £'  or  if  it'  can  be  obtained  from  jt  by  a  sequence  of  simple  transpositions. 

For  a  vector  xe  r",  we  define  xir  to  be  the  vector  (x  ^ .  ,x  ,  O  •  We 

denote  by  x  the  vector  obtained  from  x  by  arranging  the  components  of  x  in 

increasing  order.  We  say  that  x'  is  a  transposition  of  x  if  x  =  xj  and  x'  *  ct', 

where  £  ^  We  write  x  1  x".  We  note  that  this  defines  a  partial  ordering 

of  r".  This  partial  ordering  has  been  studied  by  Savage  (1957),  Lehmann  (1966), 

and  Hollander,  Proschan,  and  Sethuraman  (HPS)  (1977),  among  others. 

Let  (x,y)  e  R*' x  The  orbit  of  (x.y)  is  the  set  0.  »  {  (xir.  yo) :  it.  ocS.}. 

For  a  vector  x  €  r”  the  orbit  of  x  is  defined  similarly. 


2.1  Definition.  Let  (j,  and  (x',y')  be  two  elements  of  R.”  x  r”  belonging  to 
the  same  orbit.  We  say  that  (x,j^)  is  more  similarly  arranged  than  (x'*;^'’)  if 

^  t  a 

there  exist  Xj  2  e  such  that  xir  «  x'o  =  x  and  ;qt  ^  .  We  write  (S.;^)  2  (2',y'). 

We  refer  to  this  partial  ordering  of  1^  x  as  the  arrangement  ordering. 

&  ft  & 

We  write  (x,y)  »  if  (x,y)  ^  (x',y')  and  (x',y')  2  (x,y).  From  the  definition 

AW  ^AW  '^AW  '^AW 

it  is  clear  that  once  the  arrangement  ordering  has  been  defined  on  the  subset 
Hx,  yw) :  IT  e  S  }  of  0  it  is  completely  determined  for  all  of  0  . 

Figure  2.1  illustrates  the  arrangement  ordering  when  x  ■  (.5,  1,  3)  and 
2  >  (2,  3.5,  4).  An  arrow  in  the  diagram  from  an  element  (x,;^)  to  an  element 
(x,y')  means  that  (x,;^)  I  (x,;^') . 


((.5.  1.  3).  (2,  3.5.  4)) 


((.5.  1.  3).^ 

((.5.  1,  3).  (3.t77r^ 


(2,  4.  3.5)) 
"(r5rTf3)^(4.  2,  3.5)) 


((.5,  1.  3).  (4.  3.5,  2)) 

Figiire  2.1  .  An  Illustrative  Arrangement  Ordering. 

2.2  Remark.  Let  (x.y)  denote  the  largest  element  of  its  orbit  in  the  arrangement 

a 

ordering,  that  is,  (Ui^)  2  (xjr,;^  for  all  jr,  g  e  Then  it  is  easy  to  see 

that  -  Xj)  (y^^  -  y^)  s  0  for  all  pairs  i,  j.  In  this  case  we  say  that  x  and 
y  are  similarly  arranged.  (HLP  (1952)  use  the  expression  "similarly  ordered".) 

..  s 

We  write  x  » 

We  next  consider  the  classes  of  functions  introduced  by  HPS  (1977)  which 
are  order- preserving  with  respect  to  the  transposition  ordering  and  the  arrange¬ 
ment  ordering. 

2.3  Definition.  (i)  A  function  f  from  r”  into  R  is  said  to  be  decreasing  in 
transposition  (DT)  if  x  2  x'  ioqplies  f(x)  2  f(x')  for  all  x  e  r”.  (ii)  A  function 
f  from  R  *  R  into  R  is  said  to  be  arrangement  increasing  (AI)  if  (^c.y)  2  (x'.y') 
implies  f(jc,y)  2  f(x',y')  for  all  (g.y)  «  r"  x  r". 

Functions  which  are  AI  play  an  important  role  in  the  theory  we  develop. 

Their  properties  and  many  useful  applications  were  developed  by  HPS  (1977).  They 
used  the  name  "decreasing  in  transposition  on  R^xr”".  The  present  name  is  due 
to  MO  (1979) . 

The  domain  of  an  AI  function  is  sometimes  restricted  to  a  subset  Acr^xr’’^. 
When  this  is  done  it  is  natural  to  require  that  A  have  the  property  that  whenever 


•  -  *  •  **.  •»*  •  *.  •«  _*«*_■  *•  ’•  .*»  -*•  *•  .*»  _*•  .*•  _*• 


(x,y)  c  A,  we  have  €  A  for  all  v,  o  e  It  is  clear  from  Definition 

2.3  that  the  AI  property  of  functions  is  essentially  a  property  of  functions  on 


S  X  S  .  More  precisely,  fCx.y)  is  AI  on  A  if  and  only  if  for  each  fixed  pair 
n  n 

(x,y)  €  A,  f(^,^  is  AI  as  a  function  of  j,  £•  In  the  saute  way  the  DT 
property  is  a  property  of  functions  on  S^. 

HPS  (1977)  give  an  alternative  definition  of  an  AI  function.  Proposition 

2.4  below  states  that  their  definition  is  equivalent  to  Definition  2.3  .  The 
definition  of  HPS  is  more  useful  in  determining  whether  or  not  a  given  function 
is  AI. 

2.4  Proposition.  (MO,  1979).  A  function  f  from  x  r”  into  R  is  AI  if  and  only 

if  (i)  f(2,;^)  »  f(2E>XI^  e  R”  x  r“,  s  £  and  (ii)  f(i,)^  ^  f CX.Z')  * 

where 

2.5  Remark.  A  function  satisfying  condition  (i)  of  Proposition  2.4  is  called 

permutation  invariant  (PI) .  We  also  use  this  terminology  to  describe  a  function 
f  defined  on  r"  satisfying  f(xir)  «  f(x)  for  ir  e  S  .  Condition  (ii)  is  equivalent 
to  stating  that  for  fixed  x,  the  function  f^(;^)  =  is  DT. 

X  ~  ~ 

Mr  ^ 

HPS  (1977)  give  many  exaaqjles  of  AI  functions  including  a  number  of  well- 
known  densities  in  statistics. 

3.  Positive  Set  Functions  in  Arrangement.  Arrangement  Preserving  Kernels,  and 
Their  Preservation  Properties.  In  this  section  we  define  two  new  classes  of 
functions  and  establish  preservation  properties  of  these  functions  under  various 
statistical  and  mathematical  operations. 

Subsection  3.1  Positive  Set  Functions  in  Arrangement  and  Arrangement  Preserving 
Kernels.  Let  f  be  a  function  of  one  or  more  vector  arguments.  Define  the 
difference  operator  to  be: 


\  . - 

where  and  x^  differ  by  a  simple  transposition  of  the  i  and  the  j  components. 
We  drop  the  superscript  ij  when  it  is  understood. 

3.1  Definition.  A  function  f  from  r" x  r”  into  R  is  called  a  positive  set  function 

^ij  ^ij 

in  arrangement  (PSA)  if  x  >  x'  and  >■'  y,'  i  3  implies 


vJJ  f(x,y)  >0.  that  is.  f(x,y)  -  f(x,y')  -  f(x'.;^)  +  f(x',y')  2  0. 


This  is  the 


condition  needed  on  the  joint  density  of  (X.Y)  to  obtain  the  stochastic  version 
of  the  HLP  inequality. 

3.2  Proposition.  The  function  f  is  PSA  and  PI  iff  f  is  AI.  The  proof  is 
omitted  because  it  is  quite  simple. 

3.3  Definition.  A  function  K  from  (r” x  r”)  x  (r”  x  r”)  into  R  is  called  an 
arrangement  preserving  (AP)  kernel  if: 

(i)  K(u,  x;  V,  y)  is  PI  in  (u,  x)  and  in  (v,  y) ,  and 
(ii)  for  all  u,  v  €  r",  K(u,  x;  is  PSA  in  (x,  y) . 

We  show  later  in  this  section  that  under  mild  conditions  on  the  measure  m 
the  function  f  given  by  f(u,v)  =  whenever  g 

is  AI  and  K  is  AP.  This  preservation  property  of  AP  kernels  is  very  useful  in 
obtaining  stochastic  rearrangement  inequalities  for  random  vectors  (X(u) , 
which  depend  on  parameters  (u.v) . 

It  is  clear  from  the  definitions  that  the  properties  of  these  two  classes 
of  functions  derive,  as  in  the  case  of  AI  functions,  from  corresponding  properties 
of  functions  on  S^.  We  make  this  more  precise  in  the  next  remark. 


3.4  Remark.  (i)  A  function  f(x,^)  is  PSA  (AI)  in  (S.^) 
in  (iT,o)  for  each  fixed  x  and 


7 


(ii)  A  function  K(ii,x;  v,^)  is  AP  in  (u.x;  v,^)  iff  vt.^S)  is 

AP  in  (p,ir;  x.a)  for  each  fixed  u,  x,  v  and  y. 

In  Theorem  3.5  below  we  give  a  method  for  constructing  many  examples  of 
PSA  functions  and  AP  kernels  using  positive  set  functions.  Recall  that  a  positive 
set  function  is  a  nonnegative  function  ^  defined  on  a  subset  A  x  b  of  satisfying 

-  i^(Xj,y2)  -  i>(x2,yj)  +  <Kx2,y2)  ^0  for  all  x^  <  X2  in  A  and  <  y2  i*^  B. 

3.5  Theorem.  Let  <p  be  a  positive  set  fimction. 

(i)  If  fj  and  f2  are  DT,  then  f(x,jr)  =  *(fj(x)  ,f2(y))  is  PSA. 

(ii)  If  gj  and  g2  are  AI,  then  v.jr)  =  x)* 

The  proof  is  omitted  because  it  is  simple. 

We  conclude  this  subsection  with  some  useful  preservation  properties  of 
PSA  functions  and  AP  kernels.  The  proof  of  the  next  result  is  easy,  and  hence  is 
omitted. 

3.6  Proposition.  Let  K(u,x:  v,y)  be  AP  and  let  f()J,v)  be  constant  on  the  orbits 

°u  V  X.Z)  sKCu.x;  v,x)fCu,v)  is  AP. 

HPS  (1977)  showed  that  if  g  and  h  are  positive  AI  functions,  then 

=  8Ci»y)*'C5.)^)  is  an  AI  function.  Proposition  3.7  below  is  similar  to 
this  result. 

3.7  Proposition.  Let  g  be  an  AI  function  and  let  h  2  0  be  a  PSA  function .  Then 
f(x,y)  =  g(x,y)h(x,y)  is  a  PSA  function. 

^ij  ^ij 

Proof.  Let  X  2.  x'  and  ^  2  .  By  definition, 

f(x,y)  +  f(x',y")  =  g(x,y)h(2,y)  g(x',y')h(x',y') 

•  g(£.y)Ch(x,y)  ♦  h(x',y')3  (since  g  is  PI] 

^  g(S*^') th(x',y)  ♦  h(x,y')]  (since  g  is  AI  and  h  is  PSA] 

■  f(x',y)  f(x,y').  (since  g  is  PI]. 

Thus  f  is  PSA.  II 


s 


3.8  Corollary.  (i)  Let  and  g2  be  DT  and  let  h  be  AI.  Then  = 

8l(5)h(x.jr)g2(j;;)  is  PSA. 

(ii)  Let  gj  and  g2  be  PI  on  r”  and  let  h  be  PSA.  Then  f(x.^)  =  g]^(x)h(x,;^)g2(;t) 


is  PSA. 

Proof.  Pnrt  (i)  follows  from  Proposition  3.7  and  Proposition  3.5  (i) .  Part 
(ii)  follows  from  Proposition  3.7  and  the  fact  that  if  gj^  and  g2  are  PI  functions 
on  r",  then  g(x,j^)  =  gi(x)82^JC^  ^ 


Subsection  3.2  Preservation  Properties  of  Positive  Set  Functions  in  Arrangement 
and  Arrangement  Preserving  Kernels.  Next  we  present  the  main  preservation 
properties  of  PSA  functions  and  AP  kernels. 

Let  m  denote  a  measure  on  the  Borel  a-field  of  R.”  x  r”  such  that  for  all 
permutations  £  and  all  measurable  subsets  A^b  of  R*^  x  r”  we  have  m(A  x  B)  = 

®  AxB}.  We  call  such  a  measure  permutation  invariant.  We  also 
call  a  measure  m  on  the  Borel  a-field  of  R*^  permutation  invariant  if  for  all 
measurable  subsets  A  of  R*'  and  all  jr  e  S^^,  m(A)  =  m{xir:  xeA}. 

3.9  Theorem.  Let  g  be  PSA  and  let  K  be  AP.  Then  assuming  the  integral  exists 
finitely,  f(u,v)  =  //gCx.y)K(u,x;  v,jr)m(dx,d^)  is  PSA. 


Proof.  Let  and  A^  =  A^^.  We  have  A^Ayf(u,v)=//g(x,j^)A  A  K(u,x;  v,y)m((te,d^) . 

Partition  r”  x  r”  into  the  four  subsets  {(x,y):  x.  Sx.,  y.  Sy.},  {x,y)  : 

~  1  J  1  J 

Xi>x^,  y^sy^},  {(x,;^):  x^sx^,  and  {(x,y):  x^  >  x  ^ ,  yi>yj}-  Consider 

the  integral 


(3.1) 


I  /  t«Cx,y)A  A  K(u,x;  v,y) 

x.sx-.yi^y. 


gCx',)^)AyAyK(u.x';  v.j:) 


*  8^S»y'^Vv  v.y')}a(dx,dx) 


By  the  permutation  invariance  of  K  in  the  first  and  second  pairs  of  argviments 


we  may  rewrite  the  above  expression  as 


(3.2)  /  I  {g(x,y)A  A  K(u.x;  v.y)  -  gCx',jj)A  A  K(u,x;  v,^) 

x^xv^y 
i  j* 

-  sCx.)rn^^\U}i.x;  v.y)  +  g(x',y')A  A^K(u,x;  v.y)}m(dx.dj^) . 

We  notice  that  for  x^  =  x^  or  y^  =  y^,  the  integrand  is  zero.  Hence  the 
integral  in  (3.2)  is  equal  to  the  integral  in  (3.1). 

The  integral  in  (3.2)  is  equal  to 

X.  ^X.V.  5v  ~  ~~ 

1  J* 

Since  each  expression  in  brackets  is  nonnegative  by  hypothesis,  the  integral 
is  nonnegative.  Thus  f  is  PSA.  || 

In  Proposition  3.2  we  proved  that  a  PSA  function  which  is  permutation 
invariant  is  AI.  This  result  together  with  Theorem  3.9  yields  a  preservation 
theorem  for  AI  functions  which  we  prove  next. 

3,10  Theorem.  Let  g  be  AI  and  let  K  be  AP.  Then  assuming  the  integral  exists 

finitely,  f(u,\^  =  //g(x.;^)K(u,x;  v»;^)n(dx,d;^)  is  AI. 

Proof.  By  Theorem  3.9,  K  is  PSA.  Hence  by  Proposition  3.2  it  suffices  to 

show  that  K  is  PI.  Let  ir  e  S  and  let  be  its  inverse.  Then 

n 

fClil'Xl)  =  fIsCx,jr)K(wi,x;  ;^,y)m(dx,djr) 

=  //8(5»):)K(u,x;  v,^)m(dx,dy)  =  f(u,v) 

by  the  permutation  invariance  properties  of  g  and  k  and  a  change  of  variables.  || 

Theorem  3.10  is  the  main  result  of  this  section.  It  allows  us  to  obtain 

stochastic  versions  of  rearrangement  inequalities. 

Let  X  and  Y  be  a  pair  of  nonnegative  random  vectors  with  a  joint  density 
which  is  PSA.  For  31  ^  Xir  denote  the  vector  » •  •_•  *\^ii)^ '  shall 

see  in  the  next  section  that  the  joint  density  of  (X  ,Y  )  is  an  AP  function  of 


(w.x;  o,y) .  Let  K  in  Theorem  3.10  be  the  joint  density  of  X  and  Y.  Then  the 
theorem  states  that  if  g  is  AI  in  (x,y) ,  then  Eg(X,Y)  is  AI  in  the  indices  of  X 
and  X'  Another  way  of  stating  this  is  to  say  that  for  all  it  e 


Since  g(x,y) 


St.  tM 

S  8  C  X  j  , . 
n 

\  x.y.  is  AI, 
i*l  ^ 


•  *  »^n’  ‘  *  *^1^  * 

we  see  that  Theorem  3.10  yields  a  stochastic 


version  of  the  inequality  of  HLP.  We  will  show  in  Section  S  that  we  can  obtain 
stochastic  versions  of  many  other  deterministic  rearrangement  inequalities  for  a 
larger  class  of  random  vectors  (X<X)  which  contains  those  having  PSA  densities. 

In  that  section  we  will  also  illustrate  the  usefulness  of  stochastic  rearrangement 
inequalities  in  statistical  applications. 

As  a  consequence  of  Theorem  3.9,  we  have  proved  a  composition  theorem  for 
AP  kernels.  This  result,  presented  in  Theorem  3.1,1,  is  important  in  that  it 
provides  a  method  for  generating  new  AP  kernels  from  known  ones. 


3.11  Theorem.  Let  G  and  H  be  AP  kernels.  Then 

KCh.w;  x.z]  »  //G(u,x;  v,x)H{x,w;  ^,2)m(dx,dy)  is  AP. 

The  proof  is  straightforward  and  hence  is  omitted. 

We  next  present  another  important  preservation  result  which  enlarges  each  of 
the  classes  of  functions  we  are  studying.  It  states  that  each  of  the  classes  of 
AI  functions,  PSA  functions,  and  AP  kernels  is  closed  under  mixtures. 

3.12  Theorem.  Let  (0,F,ii)  be  a  measure  space  and  let  Fg(^  be  measurable  as  a 

function  of  6  e  0  for  each  £.  Define  FsfF^  y(d@).  Then  it  follows  that  (i)  if 

Fg  is  AI  for  each  0,  then  F  is  AI;  (ii)  if  F^  is  PSA  for  each  6,  then  F  is  PSA; 

(iii)  if  F.  is  AP  for  each  9,  then  F  is  AP. 

9 

The  proof  of  this  result  is  straightforward  and  is  omitted.  Part  (i)  of 


Theorem  3.12  is  due  to  HPS  (1977). 


-  11  - 

HPS  (1977)  prove  as  their  main  result  a  preservation  theorem  for  AI  functions. 
We  show  in  Theorem  3.13  below  that  this  preservation  theorem  is  a  special  case  of 
Theorem  3.10.  Thus  the  theory  we  have  developed  extends  the  theory  of  AI  functions 
developed  by  HPS  (1977) . 

3.13  Theorem.  (HPS, 1977).  Let  g  and  h  be  AI  functions.  Then 

fC!d»X)  * 

The  proof  is  straightforward  and  is  omitted. 

In  Theorem  3.14  below  we  present  a  generalization  of  the  preservation  theorem 
of  HPS  (1977)  for  AI  functions  presented  in  Theorem  3.13. 

3.14  Theorem.  Let  g  be  PSA,  h  be  AI,  and  m  be  PI.  Then 

=  /g(2.z)h(z,y)  m(dz)  is  PSA. 

Proof.  Let  A  =  and  A  *  A^^.  By  definition, 

Vy  8(i»2)h(2.y)m(d2)  * 

/{Cg(x,z)  -  g(x',2)]h(z,;^)  ♦  Cg(x',z)  -  g(x,z)]h(z,y')}m(dz)  . 

By  partitioning  into  the  sets  {z:  z^sz^. }  and  (z:  we  may  write 

the  above  expression  as 

/  [^gCi.z)  -  gCx'.£)]th(z,;^)  -  h(z ,;^')  ]  ♦  Cg(ji,2 ')  -  g(x',z')  3Ch(z ',jr)  -  h(z',y')3ni(dz) 

z.  s  Z  . 

1  J 

/  Cg(2,z)  ♦  g(x' .2 ')  -  g(jc' ,z)  -  g(x,z')  ]Ch(£,jr)  -  h(z ,^')]m(dz; 

z.  <  z . 

3 

by  the  permutation  invariance  of  h.  The  last  integral  is  nonnegative  since  each 
expression  in  brackets  is  nonnegative  by  hypothesis. 

Thus  f  is  PSA.  II 


HPS  (1977)  applied  their  preservation  theorem  to  ranking  problems  in  non- 
parametric  statistics.  In  Section  4  we  show  how  to  obtain  some  of  their  results 
under  slightly  more  general  assumptions.  In  addition,  we  obtain  new  results 
applicable  to  other  ranking  problems. 

Theorem  3.16  below  and  its  corollaries  will  be  useful  in  these  applications 
to  ranking  problems.  In  it  we  show  how  to  obtain  new  PSA  and  AI  functions  by 
composing  these  functions  with  certain  transformations  on  r".  By  a  transformation 
we  mean  a  mapping  from  R*'  into  r”  of  the  form  t(x)  *  (tj(j^ , . . .  .t^^Cx))  ,  where 
t^  is  a  function  on  r”.  We  write  -x  for  the  vector  (-Xj,-X2,. . . ,-x^) . 

3.15  Definition.  Let  t  be  a  transformation  from  R?'  into  r”.  Then  t  is  called  a 

rank- like  (reverse  rank-like)  transformation  if  it  satisfies : 

(i)  t(SJL)  =  £(x)E  for  all  j  and 

(ii)  t(}p  and  x(-x)  are  similarly  arranged. 

One  example  of  a  rank-like  transformation  (which  motivates  the  name 
"rank-like")  is  the  rank  order  transformation,  important  in  nonparametric  statistics. 
The  rank  order  transformation  is  the  transformation  r(j^  =  (rj(x)  , . . .  ,rjj(x))  , 
where  r^(^  is  the  rank  of  x^  among  Xj^,...,x^;  in  the  case  of  tied  x's,  average 
ranks  are  \ised.  The  vector  r(x)  is  called  the  rank  order  of  x.  Another  example 
is  the  transformation  f(x)  *  (f(Xj) , . . .  ,f(Xjj))  ,  where  f  is  an  increasing  function. 

If  f  is  decreasing  then  f(2)  is  a  reverse  rank- like  transformation. 

We  now  obtain  some  results  concerning  rank-like  transformations. 

3.16  Theorem.  Let  f  be  a  PSA  (AI)  function.  Let  t  and  s  be  either  both  rank-like 
or  both  reverse  rank-like  transformations  on  r”.  Then  ?(x,y)  =  f(t(x),s(y)) 


is  PSA(AI). 


'■11  IJ 

Proof.  Assume  that  f  is  PSA.  Let  x  s  5' and  'L  ^  X.' •  ** 

fCx.y)  +  f(x'.jr')  =  fCt(x).  s(jr))  +  f(t(x'),  s(jr'))  . 

^  fby  (i)  of  Definition  3.15] 

fCt'(x),  s(;^))  ♦  f(t(x),  s'(;|r))  »  [by  (ii)  of  Definition  3.15] 

fCt(x').  sC;^))  *  f(t(x).  s(y-))  =  f(x^y)  +  ICx,;^') . 

Thus,  f  is  PSA. 

Now  assume  that  f  is  AI.  By  Proposition  3.2  f  is  PSA.  Hence,  by  what  we 
have  shown  above,  f  is  PSA.  To  show  that  f  is  AI  it  remains  to  show,  by  Proposi¬ 
tion  3.2,  that  f  is  PI  .  Let  v  «  S  .  Then  for  all  x,  y  e  r”  we  have 

fCH.jyi)  *  fCt(xir),  s(jnr))  »  f(t(x)2,  s(jr)it)  »  f Ct (x)  ,  s(jr))  =  f(x,jr)  . 

Thus  ?  is  AI.  II 

Using  the  examples  of  ranklike  transformations  mentioned  previously,  we  can 
obtain  results  concerning  rank  orders  and  increasing  (decreasing)  transformations. 

3.17  Corollary.  Let  f  be  a  PSA(AI)  function,  (i)  If  r(x)  and  s(;^)  are  the  rank 
order  transformations  of  j  and  then  f(r(x),  s(y))  is  PSA(AI) .  (ii)  If  g  and  h  are 

either  both  increasing  or  both  decreasing  functions,  then  f(g(x),  h(y))  £ 

f((g(Xj),  g(x2),  ...,  (hCy^).  ^(yj),  h(yj^)))  is  PSA(AI)  . 

Since  the  transformation  t(3c)  =  ^  clearly  rank-like,  we  obtain  immediately 

3.18  Corollary.  Let  t  be  a  rank-like  transformation  and  let  f(jc,j^)  be  PSA(AI) . 

Then  f(x,t(y))  is  PSA(AI) . 

Another  useful  corollary  of  Theorem  3.16  is: 

3.19  Corollary.  Let  t  and  s  be  either  both  rank-like  or  both  reverse  rank-like 


4.  Stochastically  Similarly  Arranged  Pairs  of  Random  Vectors.  In  the  first  sub’ 
section  of  this  section  we  introduce  the  notion  of  stochastically  similarly  arranged 
pairs  of  random  vectors. 

We  prove  that  if  stochastically  similarly  arranged  (SSA)  then 

f(XjY)  f(X  ,YJ  for  all  A1  functions  f,  where  (X  ,Y  )  is  the  random  vector 
~  ~JL  ^  ns  "TZ 

^^ir(l)  ’  ■  ■  ■ ’\(n) ;  ’  ■  *  * ’^o(n)^  ‘  ***  ®  joint  density 

which  is  PSAorAP  then  (X.,X)  stochastically  similarly  arranged.  In  the  case 
where  X  and  X  depend  on  parameters  and  we  define  the  slightly  more  general 
notion  of  (X,X)  being  stochastically  arranged  (SA)  like  .  We  prove  results 

analogous  to  the  above  for  these  pairs  of  random  vectors. 

In  Subsection  4.2  we  show  that  many  multivariate  densities  of  interest  in 
statistical  practice  govern  pairs  of  random  vectors  which  are  SSA  or  SA  like 
parameters  (a,g)  by  showing  that  they  belong  to  the  class  of  PSA  functions  or 
AP  kernels. 

We  show  in  Subsection  4.3  that  under  certain  statistical  operations  on  pairs 
of  SSA  random  vectors,  the  property  of  being  SSA  is  preserved.  (As  an  example, 
we  show  that  this  property  is  preserved  under  mixtures.) 

We  conclude  the  section  with  some  generalizations  of  the  work  of  HPS  (1977). 


Subsection  4.1  A  Stochastic  Arrangement  Ordering  for  Pairs  of  Random  Vectors. 

We  begin  this  section  with  the  definition  of  a  new  notion  of  stochastic  comparison 
between  pairs  of  random  vectors.  This  notion  leads  to  a  stochastic  version  of  the 
arrangement  ordering  for  pairs  of  n-ttqsles  presented  in  Definition  2.1  . 

4.1  Definition.  Let  (X,Y)  and  (W,Z)  be  pairs  of  random  vectors.  We  say  that 
(X,Y)  are  stochastically  more  similarly  arranged  than  (W,Z)  if,  for  each  AI 
function  f  we  have  Ef(;!C,Y)  2  EfOJ,^),  or,  equivalently,  f(2C,X)  f(!!*S)  •  We 

write  this  as  (X,X)  (jjf,g)  . 


We  frequently  consider  pairs  of  random  vectors  (X(a),Y(B))  =  (Xj^(Oj^) , . , . , 

associated  vector  parameters  a  As  noted 

above,  we  write  (X  ,Y  )  for  (X,., , , . . .X  ,  .  ;  Y^,,. , . . . ,Y^.  J .  In  this  case  we 
can  think  of  the  indices  v  and  a  as  being  the  parameters  of  the  random  vectors 
^  and  Y  .  In  Theorem  4.2  below  we  impose  a  condition  on  the  joint  distribution 
of  X  and  Y  so  that  (a,£)  I  (a'.jl')  implies  that  (X(a),Y(jy)  *1*  (X(a')  ,Y(g'))  . 
Thus  the  arrangement  ordering  of  the  parameters  is  reflected  in  a  stochastic 
arrangement  ordering  of  the  random  vectors. 

We  use  the  notation  %  for  the  vector  of  order  statistics  (X^^ , . . . . 

4.2  Theorem.  Let  CiCH  ,X.CS)^  *  pair  of  random  vectors  with  parameters  o  and 

£.  Let 

(4.1)  P(;C(2)  =  »T,  Y(£)  »  ^llCjs)  *  X,  |()|)  .  i) 

be  a  version  of  the  conditional  probability  which  is  an  AP  function  of  (a.w;  j|,a) 
for  each  fixed  pair  (E,^)  •  T^en  (a»|)  ^  (2'»fi')  implies  (XC®) 

Proof.  The  result  is  a  direct  consequence  of  Theorem  3.10.  || 

In  the  special  case  where  we  interpret  the  indices  of  ^  and  Y  as  parameters, 
we  have  the  following  corollary  of  Theorem  4.2  . 

4.3  Corollary.  For  a  pair  of  random  vectors  (X,p ,  let 

(4.2)  P(X*|l.  X  =  ^1S  =  I» 

be  a  version  of  the  conditional  probability  which  is  a  PSA  function  of  (^r,^)  for 
each  fixed  pair  (x,^).  Then  (£,t)  I  (£',t')  implies  (X  ,Y^)  (X 

Theorem  4.2  and  Corollary  4.3  motivate  the  next  two  definitions. 

4.4  Definition.  We  say  that  (X(a) ,  Y(6))  are  stochastically  arranged  like  (a>B) 

(SA  like  (a,^))  if  the  random  vectors  (X(ig)»Y(fi))  ^  version  of  the 

conditional  probability  in  (4.1)  is  AP. 


We  will  sometimes  drop  the  parameters  and  write  (X,Y)  for  (X(a)  ,Y(3)) . 

If  F  is  a  c.d.f.  governing  a  pair  of  random  vectors  which  are  SA  like  (a,B}  we  say 
that  F  is  SA  like  (a,g). 

4.5  Definition.  We  say  that  (X,Y)  are  stochastically  similarly  arranged  CSSA)  if 
the  random  vectors  (X,Y]  are  such  that  a  version  of  the  conditional  probability  in 
(4.2)  is  PSA. 

We  say  that  the  c.d.f.  governing  a  pair  of  SSA  random  vectors  is  SSA.  It 
can  be  shown  that  if  the  pair  of  random  vectors  (^,)Q  is  degenerate  at  (£,;^ ,  then 
(}(,X)  SSA  iff  2  and  ^  are  similarly  arranged.  Hence  the  notion  of  SSA  random 
vectors  is  a  stochastic  generalization  of  the  deterministic  notion  of  similarly 
arranged  pairs  of  vectors. 

4.6  Remark.  Note  that  if  (X(op ,  Y(j|))  are  SA  like  (o,^) ,  then  for  (a,|)  fixed, 

(X(a) ,  XCf))  SSA.  Conversely,  if  (X,p  are  SSA,  then  (X^>Y^)  are  SA  like  (g,*T) . 

In  Section  S  we  use  Theorem  4.2  and  Corollary  4.3  to  obtain  numerous  stochastic 
rearrangement  inequalities  for  pairs  of  random  vectors  which  are  SA  like  (a,B)  or 
SSA.  We  show  that  these  inequalities  contain  as  special  cases  their  deterministic 
counterparts. 

In  the  next  theorem  we  present  a  sufficient  condition  for  a  pair  of  random 
vectors  (X,Y)  to  be  SSA. 

4.7  Theorem.  Let  (X,Y)  have  a  PSA  density  with  respect  to  a  permutation  invariant 
measure.  Then  (X,Y)  are  SSA. 

The  proof  is  straightforward  and  is  omitted. 

Similarly,  we  have  a  sufficient  condition  for  (X(s^,X(||))  to  be  SA  like 
(S(,J^  which  we  state  next. 

4.8  Theorem.  Let  (X(oi),Y(6))  have  an  AP  density  with  respect  to  a  permutation 
invariant  measure.  Then  (X(a),Y(B))  are  SA  like  (o,S) . 

The  proof  is  straightforward  and  hence  is  omitted. 


In  the  next  theorem  we  give  a  sufficient  condition  for  the  distribution 


function  of  a  pair  of  random  vectors  to  be  PSA. 

4.9  Theorem.  Let  fCjc,;^)  be  a  PSA  density  and  let  F(x,;^)  be  its  corresponding 
distribution  function.  Then  is  PSA. 

Proof.  Write 


F(2.X) 


// 


n 

n  I 
i»l 


(o.x^) 


(Ui) 


n 

H  If  ^ (v.)f(u,v)mCdu,dv) . 
i»l  ^  ~ 

n 


n 


Since  each  product  of  indicator  functions  is  AI,  II  I.  I^^i^  **  ®  y 

is  AP  by  Theorem  3.5  (ii) .  Since  f  is  PSA,  F  is  PSA  by  Theorem  3.9. 

An  analogous  result  holds  when  (X,Y)  has  an  AP  density. 

4.10  Theorem.  Let  f(u,x;  v,;^  be  an  AP  density  and  let  F(u,x;  v,;^  be  its  corre¬ 
sponding  distribution  function.  Then  FCju,^;  v>;^  is  AP. 

The  proof  of  Theorem  4.10  follows  from  Theorem  3.11  in  a  manner  similar  to 
that  of  Theorem  4.9,  and  hence,  is  omitted. 

The  results  in  this  subsection  allow  us  to  obtain  stochastic  rearrangement 
inequalities  for  a  large  class  of  random  vectors  which  contains  those  pairs  (jC,X) 
having  PSA  or  AP  densities.  Many  examples  of  such  densities  are  given  in  the  next 
subsection. 

Subsection  4.2  Densities  Governing  Random  Vectors  Stochastically  Arranged  like 
Parameters  Ca.fiJ •  The  purpose  of  this  subsection  is  to  show  that  many  multivariate 
densities  of  interest  in  statistical  practice  govern  pairs  of  random  vectors  which 
are  SA  like  C9,g).  In  the  next  subsection  we  show  how  operations  of  statistical 
interest  on  pairs  of  random  vectors  which  are  SA  like  (a,jg)  preserve  this  property. 
Hence,  from  the  basic  examples  given  in  this  subsection,  many  other  examples  can 
be  constructed. 

In  Theorem  4.11  below  we  show  how  to  construct  pairs  of  random  vectors  which 
are  SA  like  (a>g)  from  independent  random  vectors  with  Al  densities.  Since  a  large 
number  of  well-known  densities  are  AI,  this  result  provides  many  examples  of  pairs 
of  random  vectors  which  are  SA  like  By  Remark  4.6,  corresponding  to  each  one 

of  these  examples  we  can  obtain  an  example  of  a  pair  of  random  vectors  which  are  SSA. 


4.11  Theorem.  Let  ^  and  X  be  independent  random  vectors  in  R  x  r“  each  having 
an  AI  density.  .Then  the  joint  density  of  (X,Y)  is  AP. 

Proof.  The  result  follows  immediately  from  Theorem  3.5  (ii)  with 
♦  (x,y)  »  xy.  11 

The  following  examples  of  AI  densities  can  be  used  to  construct  AP  densities. 
(See  HPS,  1977  for  proofs). 

n  u  *i 

4.12a  Multinomial:  gj(u,x)  »  NI  y-  ,  where 

n  n 

0<Uj<l,  X.  =  0,1,2, ...  ,i  »  1,. ..  ,n,  ^  and  ^  Xj^  *  N. 

i»l  ^  i*l  ^ 

4.12b  Negative  multinomial:  g2(li»X)  * 


r  N  ♦  ^  X^ 
i«l  ^ 


?  1-N  -  [  X  "  “i  ^ 


i»l  i‘ 


where  Uj^>0,  Xj^»0,l,...,  i»l,  ....  n,  and  N  >  0 . 


4.12c  Multivariate  hypergeometric : 


u. 

n  1 
.11,  X. 

i“l  L  1 > 


^  “i 

i»l  ,  where  u.  >  0,  x.  *  0,1 , , 
N  ^  ^ 


n  n 

I  Xi*N  <  ^  u^. 

i»l  i»l 


4.12d  Dirichlet;  g^(B,2) 


0  ♦  y  u. 
i»l  ^ 


r(0)  n  r(u.) 

i*l  ^ 


n  0-1  n 

-  I 

i  i“l 


a  ,  U.-l 

0-1  n  1 
.n  *i 


where  u.  >  0,  x.  a  0,  i*l,  . . . ,  n,  ^  x.  5 1,  and  0  > 

1  X  <■  * «  X 
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4  ■  12e  Inverted  Dirichlet ;  gg  (J4,23 


®  *  “i 

i*l  J  X  i*i  1 


r(9)  n  rCu.) 

i»l 


n 

1*1 


0  n 

♦  I  *1 

i«l 


where  >  0,  i»l,  n,  and  0>O. 


4.12f  Negative  aultivariate  hypergeometric: 


N!r 


n 

^  “i, 

i-1 


n  rCx.  ♦u.) 


n  x.!r(N+  In.) 

i«l  ^  i«l  ^ 

n 

where  u.  >0,  x.  •0,1,...,N,  J  x.  ■  N,  and  . 

1  ^  i-l  *• 

4.12g  Dirichlet  compound  negative  aailtinomial;  ■ 

r(N  ♦e) 


r 

n  1 

N*  1  X. 

r 

n 

e*  i  U. 

i-1 

i-1 

n  X.  !r(N)r(9)r 

i»l  ^ 


n  n  1 
N>9^  In.*  I  X. 
i-1  ^  i-1 


n  r(x.  ♦  u.) 

*  n  - i - — 


i-1  rcu.) 


where  u^>0,  x^-0,1,...,  i-l,...,n,  9>0,  and  N*  1,2 ... .  . 
4.12h  Multivariate  logarithmic  series  distribution;  ggClJ.s) 


f  \ 

I  V' 

1  i-1  J 

f 

L-.  X 

f  ^ 

n 

i  “i 

f 

n  ] 

1-1  ‘I 

log 


1*  i  u. 


i-1 


n 

■  ^  '‘i 
i-1  ^n 


«i 


1-1  xA 


where  >  0 ,  0,.. . , 


i-1 


f  *  *  •  9 


4.12i  Multivariate  F  distribution;  gg(ll»2)  ■ 


where  0,  i »  0,1 


rCu)  n  (2u,)“^ 
i-0  ^ 


n  “i-l 

n  X. 

i«l  ^ 


n 

2  n  rcx.)(x  ♦ 
i-o  ^  ° 


n 

,n,  u»  I  X.,x.  20,  i»0,l,. 
i-1  ^  ^ 


,n  . 


4.12]  Multivariate  Pareto  distribution: 


no' 


'  " 

n 

-1 

f  «  1 

■  a(a+  1)  ...  (a  ♦  n  -  1) 

n  u. 
[i-l 

X 

y  X.  X.  -  n  1 
•  1  1  1 
i»l 

(a*  n) 


where  x^  >  >  0,  i  «  1,. . .  ,n,  and  a  >  0. 

4.12k  Multivariate  normal  distribution  with  eoiiion  variance  and  couBBon 
covariance ;  gjiCa.j)  • 

jr^e^^Cjc-ji)  J'^(jc-u)”^,  where  \  is  the  positive  definite  covariance 
matrix  with  elements  along  the  main  diagonal  and  elements  po^  elsewhere, 
p>-l/Cn-l). 

In  Theorem  4.13  below  we  give  another  method  for  constructing  AP  densities 
from  AI  densities  on  R*'  *  R® . 

4.13  Theorem.  Let  u^,  jj2.  Xj.  and  jCj  elements  of  r”  and  let  =  (aj»Sl2^ 
jc  =  (5j,223'  ^2  ***  r"*r"  and  let  g  satisfy  g(a,2)  =  gCjiE.SZ) 

for  all  V,  2  «  =  gCa.S)  Cfl(lil.Sl)  f2^a2**2^  ^ 

is  AP  in  (Uj.2j;ii2.S2^' 

Proof.  By  Theorem  3.S  (ii) ,  fj^(iij,2j^)f(a2»225  is  AP.  The  conclusion  follows 
from  Proposition  3.6.  || 

Theorem  4.13  applies  to  all  of  the  densities  listed  above  except  4.12  i,  j, 
and  k. 

Using  Theorem  4.14  below  we  can  construct  additional  examples  of  AI  densities. 
Recall  that  a  nonnegative  fmction  f(x,y)  defined  on  a  subset  A^B  of  R^  *  R^  is 


totally  positive  of  order  2  (TP^)  if  f(Xj,yj)  xf(x2,y2)  ^  f (Xi»y2)f(3C2,yj)  for 
all  *2  in  A  and  y^  ^  y2  (1977)  have  established  the  following 

relationship  between  TP^  functions  and  A1  functions. 

n 

4.14  Theorem,  (i)  Let  f(x,y)  be  TP,.  Then  g(x,y)  ■  II  f(x.,y.)  is  AI. 

^  i-1  ^  ^  n 

(ii)  Let  f(x,y)  be  a  positive  set  function  function.  Then  g(x,y)  *  \  f(x. ,y. ) 

i«l 

is  AI . 

Thus  we  may  construct  AP  densities  by  taking  one  (or  both)  of  the  functions 
g^,  i« 1,  2,  in  Theorem  4.11  to  be  the  product  of  TP2  densities.  This  corresponds 
to  the  case  where  the  conqsonents  of  one  or  both  of  the  vectors  (X,Y)  are  mutually 
independent . 

Exaiiqples  of  PSA  and  AP  densities  can  be  generated  using  densities  which  are 
Poflya  frequency  functions  of  order  2  (PF2) .  A  nonnegative  function  f (x)  is  PF2 
if  log  f(x)  is  concave  on  (>•,<•) .  Equivalently,  a  nonnegative  function  is  PF2  if 
f(x-6)  is  TP2  for  X,  0  real.  Hence  the  shifted  counterpart  of  a  PF2  density  is  TP2. 

In  Remark  4. IS  below  we  see  how  a  result  of  Karlin  (1968)  can  be  used  to 

construct  exaaiples  of  pairs  of  random  variables  (X,Y)  having  TP2  densities. 

4.15  Remark.  Often  in  reliability  situations  we  are  interested  in  components 
whose  lifelengths  have  TP2  densities  f(9,x)-  for  exaa^le,  Poisson  or  e:q}onential , 
where  9  is  a  random  variable  which  depends  on  the  environment.  Let  X  and  Y  be  the 
lifelengths  of  two  independent  coiqjonents  with  TP2  densities  F(9,x)  and  g(9,y), 
respectively.  Let  9  have  distribution  ir(9).  Then  the  joint  distribution  of  (X,Y) 
is  given  by  K(x,y)  ■  /f(9,x)g(9,y)dir(9) .  By  the  composition  theorem  for  TP2 
functions  (see  Karlin,  1968,  p.l6),  k(x,y)  is  TP2  in  x  and  y. 

4.16  Theorem.  Let  (Xj,Yj) , . . . ,  (Xj^,Yj^)  be  independent,  each  having  a  common  TP2 
density  in  (x,y)  .  Then  (X,X)  s  ((X^ , . . .  ,Xjj)  ,  (Y^ , . . .  .Y^^) )  are  SSA. 

Proof.  By  Theorem  4.14,  (j^,X)  has  a  joint  density  which  is  AI  and  hence,  PSA. 

By  Theorem  4.7,  (X,X)  are  SSA.  || 


Examples  of  PSA  and  AP  densities  can  also  be  constructed  from  Schur 


functions. 


Recall  that  a  vector  jce  r”  is  said  to  majorize  a  vector  r"  if 


^  N  11  li 

J/ci]  . ‘  “■*  J/i  •••  ^  ‘[nr 

A  function  f  is  said  to  be  Schur- concave  (-convex)  if  f(jc)  s  (2)  f(y)  whenever  x 
majorizes  y.  The  following  theorem  yields  examples  of  AI  functions  constructed 
from  Schur  functions. 

4.17  Theorem.  (HPS,  1977)  (i)  Let  f(x.j^)  -h(x-;^).  Then  f  is  AI  if  and  only  if 

h  is  Schur- concave. 

(ii)  Let  f(2,)^  =h(x+;^).  Then  f  is  AI  if  and  only  if  h  is  Schur-convex. 
Thus,  for  Theorem  4.11,  we  can  use  Schur-concave  densities  to  construct 
examples  of  PSA  and  AP  densities. 


In  the  next  subsection  we  will  study  models  in  which  observations  of  certain 
SSA  random  vectors  (X,)^  are  subject  to  contamination  (measurement  errors).  We 
will  show  that  if  the  measurement  errors  W  and  Z,  say,  have  Schur-concave  densities 
then  (}l-^W,  X'*'2)  i’S  SSA  under  some  independence  assumptions. 

MO  (1974)  give  a  number  of  interesting  examples  of  Schur-concave  densities. 

Joint  densities  of  the  form  f(s)  >  sC^^x")  where  g  is  a  decreasing  function  and 

A  is  a  positive  definite  matrix  with  X, ,  «  .. .  »  X _  and  X. .  *  X  for  i *  j  are  called 

elliptically  contoured  densities.  MO  show  that  elliptically  contoured  densities 
are  Schur-concave. 


Another  useful  example  of  a  Schur-concave  density  due  to  MO  is: 

Let  Xj,  ...,  X^  be  independent  with  PF2  density  f.  Then  the  joint  density 
n 

of  jC,  fy  ■  n  f(X. )  is  Schur-concave. 

A  •  •  X 

-V  !■! 

Subsection  4.3  Operations  On  Random  Vectors  Which  are  Stochastically  Arranged 
Like  Parameters  (a.B).  In  this  subsection  we  show  that  certain  operations  on  pairs 
of  random  vectors  which  are  SA  like  (jj.g)  preserve  the  SA  property.  First  we  show 
that  if  (X,Y)  are  SA  like  (a,^) ,  then  so  are  (r,(^,g(Y)),  where  r  and  s  are 


rank-like  functions.  In  particular  we  establish  this  result  for  the  rank  order 
of  random  vectors  which  are  SA  like  (0,8) .  Next  we  show  that  the  SA  property 
is  preserved  under  certain  contamination  models  and  give  examples.  We  also  point 
out  that  many  of  the  results  of  Section  3  can  be  used  to  perform  operations  on 
random  vectors  which  are  SA  like  (a,£)  which  preserve  the  SA  property. 

In  Theorem  4.18  below  we  show  that  if  the  random  vectors  (X(a),Y(jg))  are  SA 
like  (ot,J&)  and  r(x)  and  §_(y)  are  rank-like  transformations,  then  the  transformed 
random  vectors  (r(X(a)),  s(Y(j|)))  are  SA  like  (a,^  . 


4.18  Theorem.  Let  (XCjg)  ,  Y(g))  be  SA  like  (ci.jg) .  Let  r  and  s  be  either  both  rank¬ 
like  or  both  reverse  rank-like  transformations.  Then  (RCci)  ,S(j|))  =  (r(X(a)), 
s(Y(^))  are  SA  like  (a,g;  . 

To  prove  Theorem  4.18  we  need  a  lemma. 

4.19  Lemma.  Let  (X(a) ,  XCJ^)  be  SA  like  (a.fi)  and  let  (g(j5),§(fi))  2  (£(Ksi)). 

where  x  ®nd  §  are  either  both  rank-like  or  both  reverse  rank-like 
transformations.  Then 

(4.3)  P(R(<5)  -  u,  SCa)  *  v||(s)  =  X.  teg)  .  2) 

is  AP  in  (a,u;  jj,v)  for  each  fixed  pair  (x,2) . 


Proof.  Let  r(x)  and  s(x)  be  rank- like  transformations  of  x  and  For 
Hi  2  ^  we  have: 

PCfiCa)  S(fi)  =xll(s)  =i»  X(g)  = 

^  ^  “ij}^{s(yo)  =^ll(a)  XCg)  =x)- 

Now  the  function  Ir  ,  is  clearly  AI  in  (ji,x)  .  Hence  by  Theorem  3.S  (ii) , 
in  (u,r(x);  X.sCi;))-  By  Remark  3.4,  ^  ^  ^  is  AP  in  Cu.Ii  X.s)  ■ 


Since  (X,Y)  is  SA  like  (a,,g) ,  it  follows  from  Theorem  3.11  that  (4.3)  is  AP 


The  proof  for  the  case  where  r  and  s  are  both  reverse  rank- like  is 
similar.  || 

Proof  of  4.18  Theorem.  We  have: 

PCECa)  “an.  SCjy  =$»  ICfi)  *v)  = 

//PCRCs)  =  Ujr .  S(S)  -  w l|(a)  »  x.  ?(£)  »  J)  x  p(| c  dx,  ^ e  d£|$(a)  *  u,  |(£)  *  v) . 

By  Lenmta  4.19,  the  integrand  is  AP.  Hence  by  Theorem  3.12  (iii) , 

PCECa)  “HI.  SCft)  “  X  alica)  “  a.  1(E)  •  x)  is  AP  in  (a,ir;  g.si) 
since  it  is  a  mixture  of  AP  functions.  || 

The  techniques  used  to  prove  Theorem  4.18  and  Lemma  4.19  yield  analogous 
results  for  random  vectors  which  are  SSA.  In  Theorem  4.20  we  state  the  counter¬ 
part  of  Theorem  4.18  for  SSA  random  vectors. 

4.20  Theorem.  Let  (j^,X)be  SSA  and  let  x  snd  g  be  either  both  rank- like  or  both 
reverse  rank-like  transformations.  Then  (jg,§)  =  (jf(E)»  £(X))  SSA. 

For  random  variables  denote  by  the  rank  of  among  X^,...,X^. 

The  random  vector  R  s  (R,  ,...,R  )  is  called  the  rank  order  of  (X, ,...,X.).  Theorems 
^  I  n  ■  '  1  n 

4.18  and  4.20  have  the  following  ii^portant  corollary 

4.21  Corollary.  Let  (X,X)  be  SSA  (SA  like  (g,^).  Let  g  be  the  rank  order  of  X 
and  let  S  be  the  rank  order  of  Y.  Then  the  random  vectors  (R,S)  are  SSA  (SA  like 

(2.J8))- 

Next  we  prove  a  result  which  is  useful  in  the  study  of  contamination 
(measurement  error)  models.  In  this  type  of  model  we  are  interested  in  a  random 
vector  but  are  able  to  observe  only  a  vector  where  K  is  a  vector  repre¬ 

senting  measurement  error. 

4.22  Theorem.  Let  (X,Y)  have  a  density  which  is  PSA  (AP  with  parameters  (o,(jy). 

Let  jV,  2  be  mutxially  independent  and  independent  of  (X.X) .  ®®®b  having  a  Schur- 
concave  density.  Then  (X+R,  X*£)  SSA  (SA  like  (gijy). 


Proof.  We  prove  the  theorem  in  the  case  where  (3C,X)  has  a  PSA  density  g. 

Denote  the  densities  of  W  and  Z  by  h^^  and  h-,  respectively.  Then  the  density 

A  A 

of  (M.  is: 

fCSl.i)  *  /g(2.j:3h^(j(i-iS)h2(£-y)m(dx,d;jr) . 

Since  h^^and  h^  are  Schur- concave,  by  Theorem  4.17(i),  hj^Cw-^)  and  h2  U-X) 

MM  MW 

are  AI  functions.  The  conclusion  follows  from  Theorems  3.9  and  4.7. 

In  the  case  where  g  is  an  AP  density,  the  result  follows  in  a  similar  way 
from  Theorem  3.11  and  Theorem  4.8  .  || 

In  a  contamination  model  of  the  type  mentioned  in  the  coaments  preceeding 
Theorem  4.22,  it  is  often  assumed  that  the  measurement  error  ]jF  is  exchangeable 
multivariate  normal  with  mean  zero  or  multinomial  with  equal  cell  probabilities, 
or,  possibly  a  vector  of  independent  random  variables  with  coimnon  Poisson  distri¬ 
bution.  In  Corollary  4.24  below  ]V  and  ^  may  have  any  of  these  distributions  or 

any  of  the  ones  given  in  Exanqile  4.19  of  D'Abadie  (1981). 

To  prove  Corollary  4.24  we  need  the  following  counterpart  of  Theorem  3.12  for 

distributions  which  are  SSA  (SA  like  (g,g)). 

4.23  Theorem.  Let  (0,f,ii)  be  a  measure  space  and  for  each  0e  6  let  be  SSA 

(SA  like  (<g^,jg)).  Assume  that  for  each  x,;^,  Fg(x,;^)  is  measurable  as  a  function  of 
8.  Then  F(ji,j^)  = /Fg(2,jr)u(de)  is  SSA  (SA  like  (o,8))  . 

4.24  Corollary.  Let  (X,Y)  have  a  density  which  is  PSA  (AP  with  parameters  («,£)). 

Let  W(^  have  any  of  the  densities  in  Examples  4.19a  through  4.19g  of  D'Abadie  (1981) 
or  let  W^,...,W^  (Z^,...,Z^)  be  independent  Poisson  random  variables  with  coimnon 
intensity  parameter  X.  Then  (X+W»  fs  SSA  (SA  like  (a,g)). 

Proof.  Since  the  multinomial  and  the  multivariate  normal  random  vectors  in 
the  hypothesis  of  the  theorem  have  Schur-concave  densities  the  conclusion  follows 
in  the  first  case  directly  from  Theorem  4.22.  In  the  second  case,  trtiere  one  or 
both  of  W^,...,W^  and  Z^,...,Z^  are  independent  Poisson  random  variables,  the 
result  follows  from  the  first  case.  Theorem  4.23,  and  the  fact  that  independent 
Poisson  random  variables  with  common  intensity  parameter  when  conditioned  on  their 
sum  become  multinomial  with  equal  cell  probabilities.  || 
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4.27  Remark.  As  noted  at  the  end  of  Subsection  4.5,  the  product  of  the  sane 
PF^  densities  is  Schur- concave.  Hence  in  Theorem  4.22  one  (or  both)  of  the  con> 
taninating  random  vectors  can  be  a  vector  of  independent  PF2  variables. 

In  general,  many  results  in  Section  3  can  be  used  to  construct  random  vectors 
which  are  SSA  or  SA  like  (2>g).  Theorem  4.23  for  mixtures  illustrates  this.  As 
another  exanqjle.  Theorem  3.15  has  the  following  counterpart  for  distributions 
which  are  SSA  (SA  like  (a,g))  representing  a  slight  generalization  of  Theorem  4.22. 

4.26  Theorem.  Let  ()C,Y)  have  density  f  with  respect  to  a  permutation  invariant 
measure  m  given  by  f(2,x)  = »  “here  g^  and  g2  are 
AI  and  h  is  PSA.  Then  (X,Y)  are  SSA. 

4.27  Remark.  Theorem  4.26  has  a  reliability  interpretation.  Suppose  that  g^  and 
g^  are  densities  of  component  lifelengths.  We  may  think  of  the  parameters  and 

jg  as  random  vectors  which  depend  on  the  environment.  We  formalize  the  notion  that 
and  g  are  dependent  by  supposing  that  the  joint  density  of  3  and  jg  is  PSA.  Ihen 
by  the  theorem  we  have  that  ()l,X)  are  SSA. 

The  next  result  gives  another  sufficient  condition  for  the  distribution 
function  of  a  pair  of  random  vectors  to  be  PSA.  Recall  that  (X,Y)  is  right  comer 
set  increasing  (RCSI)  if  PCX  >x,  Y>ylx>x;Y>  y'l  is  increasing  in  x'  and  y' 
for  each  fixed  x  and  y.  Barlow  has  shown  that  this  is  equivalent  to  F(x.y)  being 
TP2.  Consequently,  the  following  theorem  is  true. 

4.28  Theorem.  Let  (X^.Yj^)  , . . .  ,(Xj^,Yj^)  be  RCSI.  Then  the  joint  siirvival  function 

of  and  X.  PU.X)  »  5  F,(x.,y.),  is  PSA. 

i.l  ^  ^  ^ 

Subsection  4.4  A  Generalization  of  a  Theorem  of  HPS.  We  conclude  this  section 
with  some  generalizations  of  the  wozic  of  HPS  in  the  spirit  of  the  previous  sub¬ 
sections  of  this  section. 

HPS(1977)  consider  a  random  vector  XCsl)  indexed  by  a  parameter  which  has 
the  property  that  its  density  f(i,£)  with  respect  to  a  pezmutation  invariant 
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measure  is  an  AI  function.  They  give  many  examples  of  random  vectors  of  interest 
in  statistics  having  the  above  property,  and  show  a  number  of  interesting  results 
in  nonparametric  statistics  for  this  type  of  random  vector.  Actually,  most  of  the 
results  in  their  paper  concerning  this  type  of  random  vector  (cf.  Theorems  4.1, 

4.4,  and  their  corollaries)  are  true  under  the  weaker  hypothesis  that 

(4.4)  P  (X  «  xirlt  3  x)  is  AI  in  a  and  ir  for  all  x. 

(That  the  latter  condition  is  weaker  than  the  requirement  that  the  density 
f(a,jt)  be  AI  follows  from  the  fact  that  P_  CX  »  *  25)  “  —  ,  whenever 

the  density  f  exists.)  SL 

We  will  say  that  the  random  vector  ^(((j)  is  arrangement  increasing  (AI)  whenever 
(4.4)  is  satisfied.  We  prove  the  following  theorem  which  yields  as  a  corollary  a 
generalization  of  Theorem  4.1  of  HPS  given  below  in  Corollary  4.32  and  a  useful 
result  concerning  AI  random  vectors. 

4.29  Theorem.  Let  ^(g)  be  AI  and  let  R(q)  s  T(XCj^) ,  where  r  is  a  rank>like 
function.  Then  ^(jj)  is  AI. 

Proof.  By  Remark  3.4  and  Theorem  3.13, 

PCBCa)  “jfflllCsO  *2)  -  I  ^Tr(xv)  “islSCsi)  “i)  is  AI  in  o  and  u. 

Now 

^o.u^2)  “  PCRC®)  *  5ll|(«)  -  Ji)  •  /PCRCo)  “ilLlit*)  *  i)  *  «  djlBCsJ  =  i)  • 

By  the  above  the  integrand  is  AI.  Since  f  (£)  is  a  mixture  of  AI  functions, 
the  conclusion  follows  from  Theorem  3.13.  || 

The  corollary  below  is  easily  seen  to  be  an  immediate  consequence  of  Theorem 

4.29. 

4.30  Corollary.  Let  X(jj)  be  AI  and  let  R(2)  be  the  rank  order  of  X(ji) .  Then 
B(S5)  is  AI. 

HPS  proved  this  result  in  the  case  where  ^((g)  has  a  density  f(2>2)  which 
is  AI. 


S.  Applications  to  Statistics  of  Random  Vectors  Stochastically  Arranged  Like 
Parameters  (a,g) .  The  theory  developed  in  Sections  3  and  4  has  ^plications  in 
a  number  of  areas  in  statistics.  In  this  section  we  present  sample  applications. 


Subsection  5.1  Stochastic  Rearrangement  Inequalities.  Using  our  theory  we  obtain 
stochastic  rearrangement  inequalities  involving  the  rearrangement  of  components  of 
random  vectors.  We  show  that  these  inequalities  contain  well-known  deterministic 
rearrangement  inequalities  as  special  cases. 

The  first  inequality  we  obtain  is  a  stochastic  version  of  the  rearrangement 
inequality  of  HLP. 

5.1  Theorem.  Let  (X,X)  be  SSA.  Then  for  all  permutations 

J,Vi  • 


i-1 


ial 


i*l 


Proof.  The  result  follows  from  Theorem  4.2  and  the  fact  that  the  function 
f(x,y)  s  xy  is  AI.|1 

Theorem  5.1  applies  to  a  large  number  of  pairs  of  random  vectors  (X,Y) .  As 
an  illixstration  let  X  =  (Xj,X2,.  ...Xj^)  and  Y  =  (Yj,Y2,. . .  ,Y^)  be  independently 
distributed  exponential,  Poisson,  or  normal  (variance  1)  with  parameters 
Oj  s  02  s  ^  Oj^  and  s  82  ^  ^  Th®"  CX,Y)  has  an  AP  density  by 

Theorem  4.14  and  hence  the  stochastic  inequalities  in  (5.1)  hold  for  (X,Y) .  Of 
course,  the  inequalities  in  (5.1)  are  true  for  any  of  the  pairs  of  random  vectors 
()C,Y)  given  in  Section  4  which  have  PSA  or  AP  densities. 

We  give  an  example  of  how  Theorem  5.1  may  be  applied  in  reliability  theory. 

This  generalizes  a  result  of  Oerman,  Lieberman,  and  Ross  (1972)  in  the  case  where 
two  vectors  are  involved. 

5.2  Example.  Suppose  that  we  have  two  stockpiles  of  n  units  each,  of  two  different 


types  of  components.  From  these  stockpiles  we  are  to  construct  n  systems,  each 
composed  of  a  component  of  type  1  and  a  component  of  type  2  arranged  in  series. 


n.  We 


A  coiq>onent  i  of  type  j  has  a  random  reliability  j»l,  2,  i»l,  ...» 
Ill  222 

assume  that  P  =  (P^>...tPjj)  and  £  =  (Pj^,  —  ,Pj^)  are  independent,  having  AI 

densities  with  parameters  i  ...  s  and  8j  s  ...  ^8^,  respectively.  Then, 

1  2 

as  we  have  seen  in  Section  4,  (£  ,P  )  are  SA  like  (a,£).  For  the  assenily  which 

pairs  the  i  component  of  type  1  with  the  ff(i)^  component  of  type  2,  the  average 

reliability  of  n  systems  is  1  r  -1 -2  .  Thus  inequality  (5.1)  states  that 

n 

i«l 

the  optimal  assembly,  in  terms  of  average  reliability  of  the  n  systems,  is  achieved 
when  the  i^^  component  of  type  1  is  paired  with  the  i^*'  component  of  type  2 .  || 


Inequalities  for  functions  of  min  (x,y)  have  been  obtained  by  Jurkat  and 


Ryser  (1966) .  They  show  that  for  nonnegative  n-tuples  x  and 


n  n  n 

n  min(x.,y.)  a  nmin(x.,y  ,.0  a  nmin(x.,y  ,  ,) 

ial  ^  ^  i-1  ^  i-1  ^ 


for  all  £  €  S„,  and 

^  U 


Mine  (1971)  obtained  the  following  similar  rearrangement  inequalities  for 
products  and  sums  of  max  (x,y) : 

n  n  n 

nmax(x,,y.)  s  n  max(x.  ,y  ,.. )  s  IImax(x.,y  .  .), 
i.l  ^  ^  i-1  ^  i-1  ^ 

and 


n 

I 

i-1 


»  n  •> 

I  max(x.,yp  s  J^»»*^*i'yv(i)5  ^  * 


n 

I 

i»l 


By  arguments  similar  to  those  of  Theorem  5.1  we  can  obtain  stochastic  versions  of 
each  of  these  inequalities: 

5.3  Theorem.  Let  (X,Y)  be  SSA.  Then  the  following  stochastic  inequalities  hold 
for  all  permutations  v: 


(5.2) 


(5.3) 


(5.4) 


St.  n 


St.  n 


n  min(X  ,Y.)  2  n  min(X  .Y  ,..)  2  n  iiiin(X.  .Y_  .  , ) ; 

i=l  ^  ^  i-1  ^  i»l  ^ 

;  .in(x.,Yj)  n-  jjrin(x..y  )  n-  j..i^.CX^,V„.^.^) ; 


”  St .  ’f'  St , 

nmax(X.  ,Y.)  s  n  inax(X.  ,Y  ,.. )  s  n  max(X.  ,Y  .  ,) ; 

.  1  11  .  ,  1  5TUJ  1  n-i+i 


(5.5) 


St.  n 


St.  n 


I  I  max(X..Y.)  s  I  max(X  .Y^^^)  s  J  »ax(X  Y  ) . 

1*1  1*1  ^  1*1 

Mine  (1971)  has  shown  that  if  x  and  y  are  nonnegative  real  n-tuples  then 

I  Cx  .y  )  s  It  (Xi.y,(j,)  X  .1 
1=1  1*1  1»1 


for  all  permutations  2  £  S^. 

n 

Since  i^(x,y)  =  -n  (x.  is  an  AI  function,  from  Tlteorem  4.2  we  obtain  the 

i«l  *  * 

following  stochastic  version  of  the  above  inequalities. 

5.4  Theorem.  Let  (X,X)  be  SSA.  Then 

J/Xj.Vj)  ’i  *<,„))  "■ 

for  all  permutations  w. 

London  (1970)  generalized  the  results  of  HLP  and  Jurkat  and  Ryser  to  obtain 
rearrangement  inequalities  for  sxims  and  products  of  functions  having  some  convex 
properties . 

Let  x>0  and  yaO.  Define  f(x,y)  =  f(l  ♦  ^x)  where  f(e*)  is  convex  for  22O 
and  f(l)  s  f(z)  for  z  2  1.  Define  g(x,y)  s  g(^x)  where  g(z)  is  convex  for  ziO  and 
g(0)  ^  g(z)  for  z2  0.  London  has  shown  that  f  and  g  are  positive  set  functions. 
Applying  Theorem  4.2  we  can  obtain  stochastic  rearrangement  inequalities  for  subs 
involving  functions  of  the  form  f  and  g. 


5.5  Theorem.  Let  (X,Y)  be  SSA.  Let  f  and  g  satisfy  the  conditions  stated  above. 


Then  n  n  Y,,..,  st  "  ^ 

sc.  p  ...  ir(ij/v  \  V 


I  f(i  ♦  Vx.)  s'  I  f(i  >  ’'^‘Vxp  s'  I  £ci  ♦  "■^■'Vx  )  . 

i“l  i»l  i=l 


and 


St*  2  ^n-i+1. 


I  gC^X  )  *s'  I  )  ’s'  I  gC  ^'^'^Vx.) 

i»l  i»l  ^  i»l  ^ 


As  an  example  of  a  ftaiction  statisfying  the  conditions  on  f,  take  f(z)  =  log(z-«-l) 
The  function  g(z)  =  z  satisfies  the  conditions  on  g.  Another  example  is  the  function 
g(z)  s  z  log(z+l). 

There  are  many  other  examples  of  deterministic  rearrangement  inequalities 
involving  AI  functions  in  the  literature.  Using  our  theory  we  can  obtain  stochastic 
versions  of  all  of  these  inequalities  in  a  unified  way. 


5.6  Remark.  We  note  that  the  stochastic  rearrangement  inequalities  we  obtain 
contain  as  special  cases  their  deterministic  counterparts.  Let  x  and  be  vectors 
of  nonnegative  numbers.  Let  A  «  :  £  e  S^}  and  B  »  {Jv  :  ;|r  e  Let  (j5,X) 

be  degenerate  at  (2>^*  l^ien  (X,Y)  has  a  PSA  density  with  respect  to  the  counting 
measure  on  A^B,  and,  clearly.  Corollary  4.3  yields  for  any  AI  function  f, 

^  ^  ^he  deterministic  rearrangment  ineqxjtality. 


Subsection  5.2  Applications  to  Tests  of  Hypothesis.  Let  (X,Y)  be  SSA  like  parameters 
•  Let  be  a  fixed  vector  qf  r"  in  the  orbit  of  g.  In  this  subsection  we 

have  developed  to  study  the  problem  of  testing  the  hypothesis 

(5.6)  Hq.*  ft  -  against  H^*’ &  ^  Oq  • 

Let  f  be  an  AI  function  and  define  the  test  T^  by 

(5.7)  T^(x,;^  ■  1  if  f(jc,x)  ^  '’oj  *  7  if  f(jS.Z)  ■  ''aj  "  °  otherwise. 


The  null  hypothesis  is  rejected  with  probability  if  (j.jr)  is  observed. 

Note  that  in  general  the  numbers  and  y(O^Y^I)  are  determined  to  give  size 
a  to  the  test. 

Let  B-  (a,j^  be  the  power  function  of  the  above  test  against  alternatives 
‘f 

(2>J^)  *  that  is,  =  ET^CXCa)  .  We  shall  need  the  following  definition 

(see  Barlow  et  al.,  1972,  Qiapter  6). 

5.7  Definition.  Let  e  r”  x  be  given.  A  test  T  has  isotonic  power 

against  alternative  (a,£)  i  (SS<j»J§q)  (with  respect  to  the  ordering  "|")  if  for  any 
(Si.Si)  “<*  (22.62)  b”«r“  such  that  (i^.Sj)  *  (*l'&i)  ‘ 

5.8  Remark.  It  is  a  consequence  of  Definition  5.7  that  any  test  T  which  is 
isotonic  with  respect  to  the  " t "  ordering  is  unbiased  for  testing 

(5.8)  (Oq.B)  »  against  (a^^.B)  <  (ao»Bo) ,  ^  * 

Note  that  by  Remark  2.2,  the  hypotheses  in  (5.8)  are  equivalent  to  those 
in  (5.6). 

It  follows  from  Theorem  3.10,  that  tests  of  the  form  given  in  (5.7)  are 
isotonic  with  respect  to  the  arrangement  ordering  and,  consequently,  by  Remark 

5.8  such  that  tests  will  be  unbiased  for  testing  H  against  H  .  We  state  this 
formally  in  the  theorem  that  follows. 

5.9  Theorem.  Let  (X,X)  be  SSA  like  (55^,)^.  Consider  testing  the  hypothesis 

Hq:  (So.jS)  -  (jSo'ficP  “a-  *  CjJo*g<P  •  * 

Let  f  be  an  AI  function  and  let  T^  be  the  test  given  in  (5.7).  Then  the  test  T^ 

a 

has  isotonic  power  against  alternatives  ^  Consequently,  a  test 

s  s 

based  on  T^  is  unbiased  for  testing  against  *  g. 

A  nximber  of  well-known  statistics  are  AI  functions  and  hence  can  be  used  to 


test  the  hypotheses  in  (5.6). 


S.IO  Example.  The  following  statistics  are  AI  functions: 

1.  Product  moment  correlation:  Pearson's  product  moment  correlation 
coefficient  is  given  by 


r  » 


I  (X.-  X.)  (y.  -  y.) 

iJ _ i _ _ 


f 

'  1.3  ^  1.3  ^ 


21 


k 


2.  Spearman's  p:  Spearman's  p  is  given  by 

I  Crj  -  Tj)  (sj  -  Sj) 

1.1 

P  a  ^  ■  ■  ■  I  ■  ■  ■  ■ 

1.3 

where  r^  is  the  rank  of  and  s^  is  the  rank  of  y^.  Spearman's  p  can  be  viewed 
as  the  sample  correlation  coefficient  coiqjuted  for  the  ranks. 

3.  Kendall's  t:  Kendall's  correlation  coefficient  t  is  given  by 


n»l  n 


f  “  I  I  C(Xj,x.)C(y. ,y.)  where  C(a,b)  - 

i»l  j-i*l  ^  ^  ^ 


1,  if  (a-b)  a  0 
.1,  if  (a-b)  <  0 


4.  A  general  correlation  coefficient  of  Daniels;  Daniels  (1948)  offers 

the  following  quantity  as  a  general  measure  of  correlation  where  x  and  jr  may  be 

either  the  observations  or  their  ranks: 

I  a(x.  ,xjb(y^,y.) 

. .  Ju2 - i - : - 


'  I  a2(x.,xj  I  b2(y.,y  ) 
i.j  ^  iJ  •’ 


7T 


where  a(x^,x^)  »  -a(x^ ,x^) ,  b(y^,y^)  -  -b(y^ ,y^)  and  a(Xj,x^),  b(y^,yj)  are  non¬ 
decreasing  with  increasing  rank  separation. 

S.  The  quadrant  test;  Lex  the  observations  or  their  ranks. 

2med  denote  the  median  of  x^,  ....  x^^  (y^,  ....  y^^)  •  The  quadrant 


u-  ia(x,)b(y,),  where 
i-1  ^  ^ 


statistic  is  the  sum 


and  b(y^)  is  similarly  defined. 

6.  Blomquist' s  qtiadrant  test  for  positive  association:  Blomquist  proposed 

the  following  test  for  positive  association: 

n 

where  a^  and  are  defined  in  (5.9),  ^2  “  and  b2  >>  1-b^. 

7.  Scores  tests  of  Bhuchongkul:  Bhuchongkul  (1964)  proposed  a  class  of  rank 

tests  based  on  statistics  of  the  form 

n 

0*1  A(r.)B(s.) 
i»l  ^  ^ 

where  A  and  B  are  nondecreasing  fvstctions. 

The  proofs  are  omitted. 

From  the  fact  that  the  product  moment  correlation  is  AI,  it  follows  that  a 

test  of  against  based  on  r  is  equivalent  to  a  test  of  against  based 

n 

on  the  statistic  I  x.y.  which  is  also  AI.  It  also  follows  that  a  test  based  on 
i»l  ^  ^ 

n 

p  is  equivalent  to  a  test  based  on  I 

Daniels  (1948)  correlation  coefficient  includes  as  special  cases,  r,  p  and  t 
by  appropriate  choices  of  a  and  b. 

The  statistic  o  in  (7)  above  reduces  to  Spearman's  p  by  taking  A(i)  *B(i)  > i. 

This  test  also  includes  the  normal  scores  test  of  Fisher-Yates  by  A(r^)(B(s^))  to 

be  the  r^-th  (s^^-th)  standard  normal  order  statistic. 

5.11  Example.  Some  further  examples  of  statistics  which  are  AI  functions  are 
given  below.  They  are  due  to  Savage  (1957)-. 


T-(l,  nj  ,  •••*  r  )  ®  J 

3  1  n  1  3 

where  d(a,b)  =1  if  a<b,  0  if  a  2  b 


4.  T^(l,  ...j  nj  Tjj  •••<  r  )  »  j  d(ii,  r.) 

n+1 

n 

5.  Tg(Ej,  . . . ,  Tj,  . . . ,  a  *^i 


6.  TgCBj,  ....  r.,  ....  r  )  ■  J  ,  m*  1,  ...,  n. 

i*l  i 

Next  we  give  some  sample  applications  of  Theorem  5.9  to  contingency  table 
analysis. 

5.12  Application.  Svqppose  that  we  ask  a  group  of  people  selected  at  random  to 
indicate  their  preference  for  one  among  k  different  objects.  If  the  group  consists 
of  males  (Group  1)  and  females  (Gro\ip  2),  say,  we  may  be  interested  in  whether  or 
not  these  two  groups  have  similarly  arranged  preferences.  Conditioning  on  n,  the 
total  number  of  persons  in  the  group,  the  numbers  of  preferences  (cell  frequencies) 


"ll’  "l2’  "ik’  "21’  "22’  •••’  "2k  jSl  J  ' 

observation  from  a  multinomial  distribution  K(£j^,Ii^;  P2>Jl2)  vith  probabilities 


2  k 

^  ^  n. .  «n  may  be  interpreted  as  an 


2  k 


Pll’  Pi2*  •••*  Pik’  P2I’  P22’  ••*’  P2k  * 

preferences  of  one  of  the  groups.  Group  1,  say,  are  known,  that  is,  the  vector 
(Pjj,  ...,  Pjjj)  is  a  known,  fixed  vector.  We  are  interested  in  testing  the  hypothesis 


“o*  '•t'li*  ^12*  •••*  ‘'ik-'  “  '■P2I’  P22'  •••*  P2k^ 

against 

^a'  ^Pll’  Pl2’  Pik^  ^  ^^21*  ^22’  ■••*  P2k^ ' 

As  we  have  seen  in  TTieorem  5.2,  K(£^.nj;£2 .52)  is  AP.  where  =  (Pij.Pj2..  •  •  .Pj,,) . 

P2=  (P2i'P22*'*'*P2k^*  Hi  “  t"n»"l2**  •• ^2  "  t"21 ’"22  ’ ' '  •"2k^  * 

5.9  a  test  of  against  is  unbiased  if  the  test  statistic  is  an  AI  function. 

k 

For  example,  the  test  which  rejects  ^  n,.n_.  <c  for  an  appropriate  c  is  unbiased 

i=l 

for  testing  H  against  H  . 

^  St 

5.15  Application.  A  similar  result  follows  if  we  condition  jointly  on  nj^,  the  total 
number  of  males,  and  on  n2,  the  total  number  of  females  (the  fixed  row  totals  case 
in  contingency  tables.)  It  is  straightforward  to  show  that  K(£j,aj;  P2»n2^  “ 
^^^1’~P^^E2»H2^  *  the  product  of  two  multinomial  densities.  As  we  have  seen  in 
Theorem  4.12,  K  is  AP.  Hence,  in  this  case,  we  also  have  that  a  test  of  against 
based  on  an  AI  test  statistic  is  unbiased. 

5.14  Remark.  In  Theorem  4.18  we  showed  that  if  (JC,X)  are  SA  like  (a,^  then  their 
rank  order  (R,^  is  SA  like  (o,£) .  Thus  Theorem  5.9  also  holds  for  test  statistics 
Tj  based  on  the  rank  order  of  (X,Y) .  A  useful  application  of  the  above  remark  arises 
in  testing  for  the  existence  of  positive  dependence  between  two  time  series.  An 
example  is  described  below. 

5.15  Application.  Studies  of  air  pollution  have  shown  that  automobile  exhaust  is 
the  major  source  of  lead  elemental  air  pollution  in  many  urban  areas.  It  is  believed 
that  automobile  exhaust  is  also  the  major  source  of  bromine  pollution  in  the  atmo¬ 
sphere.  For  a  particular  city,  we  wish  to  determine  whether  automobile  exhaust  is 
the  predominant  source  of  both  of  these  two  pollutants  or,  alternatively,  whether 
other  sources  are  responsible  for  bromine  pollution.  Suppose  that  the  concen¬ 
tration  of  lead  at  time  i,  i*l,  ...,  n,  is  known.  Let  (Xj, . . .  ,Xj^) . 


To  help  in  distinguishing  between  the  two  alternative  hypotheses,  we  test 
s  s 

H^:  ^  £  against  H^:  ^  *  J^.  where  6^^  is  the  true  concentration  of  bromine  at 

time  i,  i=  1,  ....  n,  and  6=  (6^ , . . .  ,6^^) .  Rejection  of  would  indicate  that 
sources  other  than  automobile  exhaust  contribute  to  the  bromine  pollution. 

Observations  L  on  lead  and  B  on  bromine  are  assumed  to  be  governed  by  a  joint 
AP  density  with  parameters  •  By  Theorem  5.9,  we  conclude  that  a  test  using 

an  AI  test  statistic  based  on  the  ranks  of  L  and  ^  isotonic  and  ^ ,  consequently, 
unbiased  against  H..  Nonparametric  tests  for  this  type  of  co-movement  between 
time  series  have  been  proposed  by  Moore  and  Wallis  (1943)  and  Goodman  and  Grunfeld 
(1961) . 

5.16  Remark.  Suppose  that  the  measurement  and  H  are  subject  to  errors  and  Y 
with  X  -v  MVN(0,  E(pj^))  and  Y  -v-  MVN(0,  Z(p2))  ,  where 


for  Osp  s 1.  Since  the  density  of  both  X  and  Y  is  Schur-concave  by  Corollary  4.24, 
(L+X,  B  +  Y)  is  SA  like  (X,^  and,  as  before,  a  test  using  an  AI  test  statistic 
based  on  the  ranks  of  L X  and  B  +  Y  is  isotonic. 

A  similar  result  may  be  obtained  for  the  examples  given  above  involving 
contingency  tables.  In  this  case  we  suppose  that  the  observations  and  N2,  say, 
are  subject  to  errors  X  and  Y  each  having  a  Poisson  distribution.  By  Corollary 
4.24  +  X  and  II2  X  have  AP  densities. 
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